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ABSTRACT 

The population of low-luminosity (< 10^^ erg s"^) X-Ray Binaries (XRBs) has been investigated in our 
Galaxy and M31 but not further. To address this problem, we have used data from the Chandra X-Ray Obser- 
vatory and the Hubble Space Telescope to investigate the faint population of XRBs in the grand-design spiral 
galaxy M51. A matching analysis found 25 star clusters coincident with 20 X-ray point sources within 1.5'' 
(60 pc). From X-ray and optical color-color plots we determine that this population is dominated by high-mass 
XRBs. A stacking analysis of the X-ray data at the positions of optically-identified star clusters was completed 
to probe low-luminosity X-ray sources. No cluster type had a significant detection in any X-ray energy band. 
An average globular cluster had the largest upper limit, 9.23 x 10^^ erg s"\ in the full-band (0.3 - 8 keV) while 
on average the complete sample of clusters had the lowest upper limit, 6.46 x 10^^ erg s"^ in the hard-band 
(2-8 keV). We determined average luminosities of the young and old star cluster populations and compared the 
results to those from the Milky Way. We conclude that deeper X-ray data is required to identify faint sources 
with a stacking analysis. 

Subject headings: galaxies: individual: (M51, NGC 5194) — galaxies: star clusters: general — X-rays: bina- 
ries — X-rays: galaxies 



1. INTRODUCTION 

An X-ray binary (XRB) consists of a compact object such 
as a neutron star (NS) or black hole (BH), that accretes matter 
from a donor star, which is generally a main sequence (MS) 
or red giant star. Bri ght XRBs have a verage luminosities of 
~ 1036-38 ej,g g-i ( [Heinke et all200 9 ). XRBs are divided 

Ih 



into diff'erent categories based on the mass of the donor star: 
low-mass (LMXB) or high-mass (HMXB). There are other 
divisions made, which depend on system properties, that in- 
clude X-ray bursters. X-ray pulsars, millisecond pulsars, and 
microquasars. 

In LMXB systems, the do nor stars are u sually spectral type 
G-M with masses < 1.5 M© (Hynes'2010'). As a result of both 
members being lo w-mass, L MXBs typically have periods of 
minutes to days ( ^Charles & C oe 2003 ) and transfer mass pre- 
dominantly via Roche lobe overflow. Due to the high stellar 
densities in globular clusters (GCs), LMXBs form via two and 
three-bod y encounter s at a higher rate (per unit m ass) than in 
the field ([K^|T975l |ClSk|[T975| [Fabian et al.|ri 975; Poo 



ley et al."2003'). Various studies fSarazin et al.'"2003 ; J ordan 
& al. 2004 i Kim et al. 2006 ; Sivakoff' et al. 2007 ; Kundu etaT 



2007j [Humphrey & Buote 2 008nKim et al.|2009i ) found that 
~ 20% - 70% of extragalactic LMXBs are located in GCs. 
GCs that are compact, massive, bright, a nd red (metal-rich) 
show a preference f or hosting LMXBs ( S ivakoff' et al.|200T 



Paolillo et al.|[20n 
red GCs are ^ 



Specifically, these studies indicate that 

3 times more likely to host LMXBs than blue 
GCs are. Several expla nations h ave been put forth to describe 
the metallicity eff'ect. Ivanova" (^006 ) posited a process in- 
volving magnetic braking in main- sequence stars that would 
resul t in the suppressio n of LMXB formation in metal-poor 
GCs. [Maccarone et al. ( 2004 J have suggested that irradiation- 
induced stellar winds are stronger in low-metallicity stars 
because emission line cooling is not efficient, speeding up 
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LMXB evolution. This would result in shorter LMXB life- 
times and therefore a decrease in the number of LMXBs ob- 
served in metal-poor GCs. 
In our Galaxy, Chandra has identified man v low-luminosity 



10 



30-34 



erg 



Heinke et al. 



2005 ) X-ray sources 



within GCs (|Pooley et~ar 2003 ) such as millisecond radio pul 



sars and cataclysmic variables (CVs). Millisecond radio pul- 
sars {Lx ~ (1 - 4) X 10^^ erg s"^) are believed to be the final 
stage of NS-LMXB evolution whereby the NS stops accret- 
ing appreciable amounts of material and turns on as a radio 
pulsar. Studies of millisecond radio pulsars in GCs have con- 
firmed the dynamical and metallicity eff'ects on their forma- 
tion in both the radio (|Hui et al.|2010 ) and gamma-ray ( |Abdo 
et al.|2010{ |Hui et al. 2011 ) regimes. With 22 discovered mil- 



lisecond radio pulsars in GC 47 Tuc dLori mer et al.|2003] ) and 
possibly twice as many undetected (Heinke et al. [2005 1 )7 GCs 
are a hotspot for the study of LMXB evolution. 

By contrast, HMXBs have donor stars of spectral type O 
or B with orbital per iods of days, months, or possibly years 
dLewin & van der Klis 2006). HMXBs transfer mass via 



Roc he lobe overflow and /or by Bondi-Ho yle (wind) accre- 
tion (IB ondi&Hoyle|1944, Iben et al.|1995l ). Unlike LMXBs, 
HMXBs are not found to form preferentially in star clusters. 
Instead, there is an associati on between HMXBs an d star- 
forming regions, where Shtyko vskiy & Gilfanov| ([2007|) found 
a low- significance (< 2(t) constraint on the (wider)~distribu- 
tion of HMXBs relati ve to bright H I I regio ns. A study of three 
starburst galaxies by |Kaaret et al.| ( |2004| ) suggest that bright 
HMXBs are located near young clusters (YCs) of stars. As 
a result of binary-binary or three-body interactions, coupled 
with the weak gravitational attraction to the cluster core, stars 
and compact objects can be easily ejected from YCs within a 
few Myr (Portegies Zwart et al. 1999 ; Kaaret et al. 2004 ). Dif- 
ferent groups ( Kalogera 2006; Sepinsky et al. 200 5^ have used 
the StarTrack population synthesis code (Belczynski et al. 



2008|) in combination with a prescription for young cluster 



potentials to show that ^ 70% of binaries (born within young 
cluster potentials) are ejected anywhere from 30 - 100 pc as 
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a result of supernova (natal) kicks for cluster ages of 15 Myr 
or older. The fact that few bright (~ 10^^ erg s"^) XRBs are 
discovered near clusters is in line with observations of nearby 
starburst galaxies where ejections of > 200 pc (Rangelov et al. 
|2011t|Kaaretet al. 2004 ) have been found. 

The role of dynamics in the dense environments of clus- 
ters is crucial in the evolution of binaries. For instance, not 
only do ejections illust rate the interactions occurring within 
these systems but Vanb everen et al.J p012) state that bina- 
ries may have a significant effect on the chemical evolution of 
GCs. This link outlines the importance of binary evolution in 
stellar systems. Detailed theoretical modelling of binary evo- 
lution in clusters is complex and therefore requires observa- 
tions to provide much-needed constraints. The emergence of 
a connection betwee n the study of bin ary compact objects and 
gravitational waves ( [Belczynski et al.^2Q0 2 ) is another source 
of motivation for uncovering the evolutionary processes of bi- 
nary systems. 

Before the launch of the Chandra X-Ray Observatory in 
1999, the study of X-ray sources was limited to the Local 
Group due to the low angular resolution and/or sensitivity 
of previous missions. Much progress has been made inves- 
tigating XRB populations in the last decade as a result of 
Chandra 's subarcsecond resolution. However, there are still 
many unanswered questions. Currently, mostly bright XRBs 
are observed in other galaxies; M31 is the only nearby (large) 
galaxy where Chandra can observe sources below ~ 10^^ 
erg s"^ (Z hang et al. 2011) . In general, since most elliptical 
galaxies are devoid of gas and dust (and therefore lack star 
formation), only LMXBs are formed and so it is straightfor- 
ward to identify them. The younger stellar population of spi- 
ral galaxies allows the formation of HMXBs, which have been 
identified separately from the diff'use emission in these galax- 
ies via bright optical counterparts. Conversely, LMXBs have 
been more difficult to identify in spiral galaxies (because of 
faint optical counterparts) and therefore comprise an impor- 
tant element of future research. 

M51 has been studied numerous times in X-rays by dif- 
ferent groups, with the m ost recent p opulation analysis from 
Chandra observations by [Terashima & Wilson ( 20Q4| ). They 
focused on the bright-end of the X-ray luminosity func- 
tion (XLF) of M51, investigating ultraluminous XRBs, and 
implied that these were likely stellar mass BH s accreting 



near/a bove LEdd- This work wa s later extended by [Terashima 
et al.|(2006) and|Yoshida et al.|([2010) to d etect optical c oun 



terparts and variability respectively. Kilgard et al.|([2005) fur 



ther studied MSl's XLF (minimum point source luminosity 
of 2.2 X 10^^ erg s"^) in comparison to other nearby spiral 
galaxies and found a flatter slope indicative of a larger HMXB 
population. They also produced X-ray color-color diagrams 
in order to classify point sources (there is controversy about 
the use of this classification due to the spectral state vari- 
ability of XRBs). Using these results, coupled with a study 
of variability, they concluded that a minimum of ~ 27% of 
all source s detected in a sam ple of 1 1 galaxies are accreting 
XRBs. Kilgard et al. ( 2006) performed a multi- wavelength 
analysis of M5rs X-ray sources by matching them to opti- 
cal data from the Hubble Space Telescope (HST). A total of 
1 1 sources were matched to O or B-type stars and are likely 
HMXBs. They found 9 sources in you ng star clusters that 
seem to be a mix of XRBs and SNRs. iMaddox et al.| ( |2007| ) 
was able to match X-ray sources with radio observat ions and 
later confirmed 6 of the 24 sources identified in Ha by |Kilgard 



et al. (2006) as SNRs. An absence of optical counterparts co- 
incident with M51b (also known as NGC5195, the companion 
galaxy interacting with M51) X-ray sources suggests that the 
population should be composed of LMXBs. Therefore, by 
using multi-wavelength data, approximately half of the X-ray 
sources were optically classified as HMXBs or SNRs, with 
the remainder obscured or faint XRBs. 

Besides deeper observations, there is a way we can probe 
the faint-end of the XLF and investigate low-luminosity 
XRBs. This is possible wit h a procedure known as stacking. 
By utilizing the methods of [Brandt et al.| ([2001b c) and Horn-' 



schemeier et al.| ( 2001 ) we will stack star cluster positions in 



an X-ray image of M5 1 to produce in eff'ect a deeper exposure 
of the average star cluster. By producing a final stacked im- 
age of many star clusters we hope to build up a source signal 
above the background level that would be evidence of faint 
sources. This will allow us to probe the low-luminosity end 
of the XLF and determine, on ave rage, if star clusters in M5 1 
host XRBs. Cha ndar et al.[P004| ) found a mean V-I color of 
M51 GCs identical to the peak of met al-poor G alactic GCs, 
where [Fc/HImw-gcs = -1.61 (Burgarel la et al.|2 001). There- 
fore, the largely blue metal-poor GC population of M5 1 does 
not allow us a contrast with red metal-rich GCs. However, this 
makes a comparison with GCs in elliptical galaxies straight- 
forwa rd since they have ne ar identical metallicities to MSl's 
GCs ( Chandar et al.[[2QQ4 ). By separating our clusters into 
different categories based on color (age) we investigate both 
HMXB and LMXB populations. Moreover, this analysis will 
help us better understand the metallicity (estimated from op- 
tical colors) dependencies of XRBs in both YCs and GCs. In 
particular, observations that reveal the luminosity ranges and 
color correlations are crucial in constraining compact object 
and binary models such as the S tarTrack population synthe- 
sis code (B elczynski et al.|2008 ). More specifically, what do 
the average X-ray luminosities of YCs and GCs tell us about 
XRB populations? 

We adopt a distance to M51 of 8.4 ± 0.6 Mpc determined 
by the planetary nebula luminosity function ( Feldmeier et al. 
[1997|. This agrees almost exactly with the recent determina- 
tion of 8 .4 ± 0.7 Mpc using supemovae 2005cs and 201 Idh 
( [Vinkoetal.|2011) . 

2. OBSERVATIONS 

2.1. X-ray Data 

Chandra has observed M51 three times since 2000: 2000 
June 20, 2000 June 23, and 2003 August 7, for exposure times 
of 15.06 ks, 27.15 ks, and 48.61 ks respectively. The X-ray 
point source catalogue that we use for this analysis was taken 
from Liu ( 201 1 , hereafter LI 1) and corresponds to that associ- 
ated with the 48.61 ks observation (ObsID 3932). This obser- 
vation was obtained with Chandra's Advanced CCD Imag- 
ing Spectrometer (ACTS) instrument in (TIMED) VFAINT 
mode with the galaxy located on the back-illuminated S3 chip 
(proposal 04600406, RL: Yuichi Terashima). |L11| identified 
109 X-ray point sources. Taken with the longest exposure 
time, these data give us the most complete catalogue of X-ray 
point sources in M51 to date. For reference, in the 0.3 - 8 
keV energy band the limiting luminosity in the X-ray cata- 
logue was Lx ~ 4.6 X 10^^ erg s"^ and the maximum was 
~ 5.4 X 10^^ erg s"^ while the average of all point sources 
was Ljc ~ 2.4 X 10^^ erg s"^. 

2.1.1. Data Reduction 
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The 48.61 ks data were reprocessed with the updated Chan- 
dra Interactive Analysis of Observatio ns (CIAO) tools pack- 
age version 4.3 ( [Fruscione et al.|2006| ) a nd the Chandra Cali- 
bration database (CALDB) version 4.4.6 ( [Graessle et al.|2006| ) 
to clean the image for optimal analysis of faint sources. Start- 
ing from the level- 1 events file, acis_runJiotpix was used 
to produce a bad pixel file that identified observation- specific 
bad pixels, hot pixels, bright bias pixels, and afterglow events. 
Since the observation was completed in VFAINT mode and 
we are studying sources with a small number of counts, 
acis_detect_afterglow was used to eliminate afterglows 
with very few events. Next, we ran ac is .process .events 
to update the level- 1 events file for grade and status informa- 
tion, charge transfer inefficiency (CTI), time-dependent gain, 
and pulse height. The energy-dependent subpixel event repo- 
sitioning (EDSER) algorithm parameter (under the pixel ad- 
justment option) was used to improve the astrometric accu- 
racy of our images. 

Each CIAO tool has unique processing parameters that need 
to be adjusted based on the given dataset. Specifically, for 
VFAINT data, acis_process_events was set to use a 5 x 5 
event island to identify cosmic-ray background events as op- 
posed to the usual 3x3 block. After filtering using the stan- 
dard (ASCA) grades (0, 2-4, 6), status bits (0), and the good 
time intervals, we produced images in various energy bands. 
These include fully processed soft, hard, and full-band images 
in the energy ranges of 0.3 - 2 keV, 2-8 keV, and 0.3-8 keV 
respectively. Energies above 8 keV are not used since the ef- 
fective area of the Chandra High Resolution Mirror Assembly 
steep ly decreases with ene rgy while the background increases 
(e.g., [Brandt et al.||2001a| ). For each image we computed an 
exposure map using the CIAO tools asphist, mkinstmap, 
get_sky_limits, and mkexpmap, which are included in a 
provided script called fluximage. We used monoenergies 
of ~ 0.93 keV, ~ 4.7 keV, and ~ 2.1 keV, which were the 
average energies for the soft, hard, and full band images re- 
spectively. These energies were used to calculate the instru- 
ment map as opposed to using weighted spectrum files. Using 
the CIAO tool dmimgcalc, we divided each image [counts 
pixel" ^] by its corresponding exposure map [cm^ s counts 
photon"^] to obtain final fluxed images in units of photons 
cm"^ s"^ pixel" ^ Each image has dimensions of approxi- 
mately 9.7' X 9.7'. Using the distance to M51 of 8.4 ± 0.6 
Mpc, we determine a linear scale of 40.7 pc arcsecond"^ . The 
plate scale of the Chandra images is 0.5 arcsecond pixer\ 
which corresponds to ~ 20 pc pixel" ^. 

In order to perform our stacking analysis accurately, both 
optical and X-ray images we re registered to the 2MASS point 
source catalogue ( Skrutskie et al."2006') reference frame. A 
total of 99 reference stars were found within the field of view 
of the Chandra image from the 2MASS catalogue. Using 
the program immatch in the WCSTools package, we found 6 
reference stars that matched within a tolerance of 0.5'' of X- 
ray source counterparts. This resulted in applying astrometric 
shifts of -0.210" in right ascension and -0.073" in declina- 
tio n to t he image and to the point source catalogue compiled 
bypT] 



2.2. Optical Data 

The star clus t er cata logue used in this analysis, taken from 
Hwang & Lee ( 2008| hereafter HL08), was compiled from 



15640 s respectively (proposal 10452, RL: Steven V. W. 
Beckwith). Hereafter, B, V, and I will be used to denote col- 
ors corresponding to B/7435^, and Ipum respectively. 
The plate scale of the HST ACS images is 0.05 arcseconds 
pixel" \ which corresponds to ~ 2 pc pixel" ^ The image 
has dimensions of approximately 6.8' x 10. 5'. We use all 
2224 star clusters that have been identified by |HL08| as Class 
1, meaning they are nearly circular and have no prominent 
nearby neighbours. 

As a result of the poor astrometric accuracy of the HST 
guide star catalogue, significant shifts are introduced in HST 
images. We follow the same procedure as described above for 
determining the off'set in the X-ray image. From the 2MASS 
catalogue we found a total of 98 reference stars in the HST 
image. The HST off'set was adjusted by applying astromet- 
ric shifts of -0.116 ± 0.320" in right ascension and -0.653 
± 0.441" in declination to the star cluster catalogue compiled 
by |HL08| These values were calculated by matching 21 ref- 
erence stars within a tolerance of 0.9" of image stars. 

3. DATA ANALYSIS 

3.1. X-ray Properties of Matched Star Clusters 

In order to perform our stacking analysis, we first needed to 
determine if there were any direct matches of star cluster posi- 
tions to X-ray point sources. We performed a matching anal- 
ysis with TOPCAT (Taylor 2005 ) using the centres of each of 
the 2224 identified star clusters (HL08 ) with the 109 X-ray 
point sources identified by |L1 1 We found 20 matches to 25 
star cluster positions (multiple star clusters can be associated 
with one X-ray point source) within a 1.5" radius. We used 
a slightly larger matching radius than previous studies (gen- 
erally ~ 1") such as Sivakoff' et al. (2008 ). This value was 
used because we are observing mostly young clusters, where 
ejections are more likely, as opposed to GCs. The results are 
shown in Table |3] along with an optical image of the matches 
in Figure fl] A close-up of the star cluster groups outlined 
in Table [3]is shown in Figure [2] To assess the chance co- 
incidence probability of our matches, we follow ed a method 
similar to that used by Antoniou et al.| ( j2009 ) and Zezas et al. 
( 2002| ). By creating 25 source lists from our original X-ray 
catalogue with each source offset by ± 5-10", we determined 
a false-match rate of ~ 50%. This could be expected with a 
large matching radius and sample size. Using the 1" matching 
radius our false-match rate is reduced to ^ 20%. A color-color 
plot of the X-ray point sources is shown in Fi gure [3| Over- 
laid is the X-ray color classificat ion scheme of Kilgard et al. 
( |2005| ), which was modified from |Prestwich et al.| ( |2003| ). The 



values in color space are outlined in Table [IJ XRBs have 
been combined into one group since delineating the popula- 
tion as low-mass or high-mass is difficult and controversial. 
Although this color classification is somewhat arbitrary (an 
XRB can end up in diff'erent regions of the diagram depend- 
ing on its spectral state), it is a good approximation to the 
types of sources we are observing. In addition, the uncertain- 
ties associated with X-ray colors are large because of the low 
counts per source, increasing the probability of possible mis- 
classification. Soft and hard colors (SC & HC) are defined in 
equations [T] and |2] respectively based on the soft (S, 0.3 - 1 
keV), medium (M, 1 - 2 keV), and hard (H, 2 - 8 keV) band 



data taken by the HST's Advanced Camera for Surveys (ACS) 
with the F435W, F555W, F814W, and F658N filters. The ef- 
fective exposure times were 16320 s, 8160 s, 8160 s, and 



energies consistent with LI 1 



SC = 



(M-S) 
{H + M + S) 



(1) 
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HC = 



(H-M) 



(2) 



The majority (~ 65%) of sources in Figure [3] are classified 
as XRBs. Figure [T] shows that most matches are located along 
the spiral arms ofM51 and are likely HMXBs. The exceptions 
are matches 17-20, which appear in M51b or are found at 
large radii from the nucleus of M51. Matches 18 and 19 may 



be bright LMXBs since Kil gard et al.| ( [2006| ) indicated a lack 
of optical counterparts in M51b. Match 19 is located 6'' from 
the co-ordinates of SN1945A^ but would not be associated 
with emission from the remnant. If we assume that the ejecta 
from SN1945A has a similar velocity to that of SN1987A, 
the SNR would have a radius of only ~ 0.5'' (natal kick to 
the compact object is negligible). The remaining sources are 
mainly classified as SNRs (matches 6, 10, and 12 - 14) and 
all reside within knots of the spiral arms, adding confidence 
to their classification. 





Fig. 2. — Close-up of the star clusters from Figure [TTthat were identified as 
groups in Table |3] From left to right, cluster group IDs from Table [3] are 1, 
2, and 3 respectively. All panels have the same dimensions (10'' across by 
6.1" high) and all circles are ?>" in diameter. The X in each panel marks the 
location of the X-ray point source matched to each cluster within 1.5''. 



48.61 ks Observation 
1.5" Matching Radius 
Sources: 109 
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Fig. 3. — X-ray color-color diagram of 109 X-ray point sources in M51 from 
|L11| Of these point sources, 20 were matched within 1.5" to 25 optically 
identified star clusters. The X-ray color classification scheme of Ki lgard et al?] 
{2005 1 (modified from Prestwich et al._2003j is shown with the values in color 
space outlined in Table[^ 



TABLE 1 
X-RAY Color-Color Classifications 



Classification 



Range 



Supernova Remnant 
X-ray Binary 
Background Source 
Absorbed Source 
Supersoft Source 



HC < -0.2, SC < -0.3 

-0.8 < HC < 0.8, -0.3 < SC < 0.6 

HC > -0.2, SC < -0.3 

SC > 0.6 

M=0, H=0 



Note. — Classification scheme of [Kilgard et al.|{2005) . 



Fig. 1. — Hubble Space Telescope ACS image of M51 in the F435W filter 
showing star clusters matched to X-ray point sources using TOPCAT. Table 
|3]has position and separation details including star cluster group information 
not shown here (see Figure[2|. 



3.2. Optical Properties of Matched Star Clusters 

The optical data for the 2224 star clusters in M51 (|HL08]) 
is used to make the color-color plot in Figure |4] A theoretical 
ev olutionary track from the Simple Stellar Population (SSP) 
of [Bruzual & Chariot] ( |2QQ3 :) with Z = 0.02 is overlaid after 
applying internal reddening of E{B - V) = 0. 1 . The internal 
reddening was computed using Ay^^^^ = 0.31 mag from HL08 , 

3 Obtained from the NASA/IPAC Extragalactic Database (NED) 



Foregro und reddening of E(B - V) = 0.035 from the dust 



maps of ISchlegel et al. dl998) has also been applied to the 
SSP model. The 25 labelled star clusters correspond to the 20 
X-ray point sources they were matched to from Figure [3] Star 
cluster colors were calculated by HL08 using a 6-pixel aper- 
ture in each band. The SNRs indicated are located in clusters 
that are younger than 100 million years with the exception of 
match 12. This is likely a result of the increased reddening 
near t he nucleus of M51 (E(B - V) = 0.25 - 0.3; |Lamers| 
[eraLl[2002 i Kaleida & Scowen ||20T0t where match 12 is lo- 
cated, which would force it to move into the younger region 
of color space (refer to reddening vector in Figure [4]). The 
same may be true for matches 2, 18, and 19, which are all lo- 
cated in heavily obscured regions. Therefore, most XRBs are 
classified as HMXBs since they are associated with clusters 
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younger than 100 Myr. Match 20 may be an absorbed back- 
ground source (AGN) because it is located far from the extent 
of the spiral arms and has red V -I colors in Figure [4] To de- 
termine whether our matches correspond to the brigntest star 
clusters, we created the color-magnitude diagram show in Fig- 
ure [5] Though previous studies of XRBs in elliptical galaxies 
have found th at they are preferentially found i n luminous clus- 
ters ( .Sivakoff et al.pOOTllPablillo et al.|2QTT] ), our matches do 
not seem to favour this relation as they are not found within 
a particular region in color-magnitude space. However, the 
matches located in the most luminous clusters would be more 
likely to be associated with an X-ray point source (because 
luminous clusters are relatively rare), whereas clusters on the 
outskirts may be spurious matches. Many of the faint clusters 
were those mentioned to reside in regions of high extinction 
and could therefore be highly obscured as opposed to false- 
matches. 

It is important to compare our matching results to previous 
studies to determine if any similarities exist among star clus- 
ter X-ray sources in different galaxies. [Sivakoff et al. (2007) 
found 270 of 6488 (~ 4%) GCs in 1 1 Virgo Cluster early-type 
galaxies with X-ray emission. This is in contrast to the 25 
of 2224 (~ 1%) of star clusters that show X-ray emission in 
M51. However, we used a matching rad ius that is 0.5'' larger 
than that used by Sivakoff et al. (2007). Using this smaller 
radius, we find that only ~ 0.5% of star clusters are matched 
with X-ray point sources. This much lower value could be a 
consequence of the differing XRB populations of spiral and 
elliptical galaxies. While we have identified mainly HMXBs 
in M5rs YCs, early type galaxies' XRB populations are com- 
prised ofLMXBsJn_GCs^_a^ more time for them to 



-0.5 



form ( |Bogdan & Gilfanov 2010*). Combined with the fact that 
HMXBs are short-lived, we might expect a larger fraction of 
star clusters with X-ray emission in elliptical galaxies. 

4. X-RAY STACKING 

To further investigate the X-ray source population, specifi- 
cally that of faint or quiescent XRBs, we performed a stacking 
analysis of the 2199 unmatched star clusters i n M51. Follow- 
ing th e methods of Brandt et al. ( 2001b c ) and |Hornschemeier 
|et al.| ([2001 ) we stacked 21 pixel by 21 pixel regions cen- 
tred on star cluster positions in X-ray images of M5 1 . This 
was completed using code written in the Interactive Data Lan- 
guage (IDL). 

4.1. Method Results 

Before proceeding with stacking, all star cluster positions 
that were closer than 1 1 pixels to the edge of the X-ray image 
were excluded because a stacked image (21 by 21 pixels) re- 
quires this minimum number in each direction from the centre 
to be complete. This left a total of 2187 star clusters in the en- 
tire sample. We then subdivided star clusters into 5 different 
groups based onB-V color as shown in Table [2] Each group 
was stacked in the full (0.3 - 8 keV), hard (2-8 keV), and 
soft (0.3-2 keV) energy bands. Our initial method consists of 
simply stacking all star cluster groups with no modifications 
to the star cluster regions in the X-ray image. The results are 
shown in Figure [6] 

The net counts m the source aperture (C„) were determined 
through the relation Cs - C^x As/Aiy, where C and A repre- 
sent the counts and area respectively for the source (s) and 
background (b) regions. Our unce rtainties were calculated 
using poisson statistics (method of |Gehrels||1986| ). We first 
computed the 3(T-clipped mean jj (Ch/Ah) in the background. 
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Fig. 4. — Optical color-color diagram of the 2224 star clusters in M51 taken 
from the catalogue of HL08 A theoretical evolutionary track from the Simple 
Stellar Population (SSP) model of Bruzual & Charlotj 2003 > with Z = 0.02 is 
included with internal reddening of E{B - V) = U.l and foreground reddening 
E(B - V) = 0.035 applied. Labels from 6.0 to 10.3 represent log(age) in 
years. Numbers from 1-20 are color-coded by X-ray type as indicated in 
the l egend and correspond to the 25 matches to 20 X-ray point sources from 
|L11| Numbers with an asterisk correspond to matches with separations < V 
irom Table |3] The foreground reddening vector is of length Ay = I mag. 
We have excluded clusters in the figure with (B - V) > 1.5 appearing well 
below the billion-year mark for illustration purposes. Matches 2, 12, 18, and 
19 may all be heavily reddened (see Figure^ and therefore we may expect 
them to move into the young region of color space. Match 20 may be an 
absorbed background source (AGN) because it is located far from the extent 
of the spiral arms (Figure|4j and is quite red in y - /. 



TABLE 2 
Stacked Star Cluster Groups 



Group Type 


Number 


Optical color Range 


All Clusters 


2187 


-0.23 <(B- 


V) < 1.63 


Very Young Clusters 


305 


-0.23 <(B- 


V) < 0.05 


Young Clusters 


1948 


-0.23 <(B- 


V) < 0.25 


Intermediate Clusters 


574 


0.30 <(B- 


V) < 0.75 


Old Clusters 


65 


0.75 <(B- 


V) < 1.63 



where any pixels with counts greater than yu + Scr were ex- 
cluded. The upper and lower limits for the source and back- 
ground regi on counts wer e calculated using equations (9) and 
(14) from Gehrelsl ( [T986 ). We used a value of ^ = 1.282 
(number of Gaussian cr) with confidence level parameters of 
JS = 0.01 and y = -4.0. The uncertainties in C„ were de- 
termined by normalizing the background values to that of the 
source aperture. Luminosities were computed from fluxes as- 
suming a distance of 8.4 Mpc. The uncertainties in luminos- 
ity were determined by taking a ratio of the luminosity in the 
source aperture to the counts in the source aperture and mul- 
tiplying that ratio by the uncertainties in counts. 
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Fig. 5. — Color-magnitude diagram of 2224 star clusters in M51 taken from 
the catalogue of HL08 Numbers from 1-20 are color-coded by X-ray type 
as indicated in the legend and correspond to the 25 matches to 20 X-ray point 
sources from LI 1 The foreground reddening vector is of length Ay = I mag. 
Numbers with an asterisk correspond to matches with separations < V from 
TableE 



Our stacking yielded non-detections for all cluster types 
and energy bands. In most cases, except the soft and full- 
band old clusters, the background- subtracted luminosities in 
the source apertures were negative. This means that the back- 
ground level is actually higher than that of the source. Inspec- 
tion of the images shows that many contain point sources dis- 
tributed randomly outside the source aperture. This could be 
consistent with ejections of XRBs from star clusters in M51. 
From the diameter of our image (~ 425 pc) it is possible that 
the point sources we see have been ejected from clusters up 
to ~ 200 pc, similar to results from other groups. This could 
contribute to the lack of both bright and faint XRBs within 
1.5'' (~ 60 pc) of clusters. 

In order to be certain that our results were not biased by 
the numerous bright point sources that are found in the back- 
ground region, we removed the image(s) from the final stack 
where these point sources were found. This was accomplished 
by producing a histogram that represented the mean number 
of counts per pixel in the entire 21x21 region for each star 
cluster stacked. In computing the mean value of each stacked 
star cluster image, the pixels in the source aperture were all 
set to zero. This allowed us to retrieve the value of the back- 
ground without excluding any star clusters that had significant 
source aperture counts. This may be somewhat skewed since 
a stacked image with a large mean value might also have a 
strong source signal. A sample histogram is shown in Figure 
|7]for young clusters. The tail in the histogram has been re- 
moved (black portion) since it constitutes the star cluster po- 
sitions stacked that have large background values and likely 
contain point sources. 



After evaluating every stacked image using this method, we 
compared the results to those in Figure [6] These new images 
looked uniform and resembled background noise. In contrast 
to the clusters from Figure [6| this approach has resulted in 
positive source luminosities for most clusters in the full and 
soft bands. However, the opposite is true for old clusters 
stacked using this new method, likely due to the eff'ect men- 
tioned above where more source signal was removed com- 
pared to the background. Nevertheless, this procedure did not 
result in detections within any of the source apertures. 

The fundamental problem with stacking outlined earlier is 
that diff'use emission in spiral galaxies makes analysis com- 
plicated. This X-ray emission comes mainly from HII regions 
and is a result of supemovae and OB and/or Wolf-Rayet stars 
heating gas in the interstellar medium to millions of degrees 
through shocks. In an attempt to reduce the diff'use emis- 
sion that may be contributing to the background level, we 
have identified and removed the 150 brightest (Ha lumin qsitv 
> 10^^ erg s ~^) HII regions taken from the catalogue of 



Lee 



et al. ( 2011 ). We do this by excluding the star clusters asso- 
ciated witn these identified HII regions from the final stack. 
Again, it is possible that by excluding a star cluster this way 
we may be eliminating signal from our source aperture in the 
proces s. W e have also applied the HST off'sets mentioned in 
section 12^2] to this HII region catalogue since it uses the same 



data. 

All stacked images produced using this method were again 
consistent with non-detections. Comparing to Figure |6] we 
found that numerous point sources and diff'use areas of X- 
ray emission were removed. This method, however, does 
not eliminate all point sources from the background as did 
the previous approach. For both all and young cluster types 
our source luminosities are now positive, which would be ex- 
pected since young clusters are generally located near HII 
regions. Conversely, we find no change to intermediate- age 
clusters. The old cluster sample soft and full-band source lu- 
minosities changed from positive to negative. Lastly, we ap- 
plied both methods mentioned above to the original stacked 
images to see what eff'ect, if any, this would have on produc- 
ing a source signal. No significant changes were observed as 
we retrieved upper limits for all stacked images. 

All the methods used above in stacking star clusters over- 
look the fact that only a small percentage of clusters may have 
faint XRBs. This would systematically average out any sig- 
nal that could be detected as star clusters with no XRBs are 
stacked. In an eff'ort to address this issue, we have used the 
findings mentioned in section[T]that state bright, compact, and 
massive clusters are known to be more likely to host XRBs. 
Using V magnitudes and eff'ective radii of M5 1 clusters from 



HL08] we stacked the brightest and most compact (assumed 
small Reff) 10% of clusters. This was done for the complete 
cluster sample as well as the various age ranges (Table[2]). The 
results were unchanged as non-detections followed for all at- 
tempts. Another issue that arises when stacking star clusters 
at large off'-axis angles is the increased size of the PSF, which 
would ultimately skew the flux calculations based on fixed 
aperture size and background region used throughout. Ideally, 
we would like to choose the aperture size and background re- 
gion based on the PSF size at the location of each star cluster 
that was used in a stacked image. However, because we com- 
pute luminosities on the final stacked images, we require all 
the apertures and background regions to be of identical size. 
We addressed this problem by stacking cluster types in vari- 
ous oflf-axis angle ranges where the PSF-size is approximately 
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Fig. 6. — Stacked Chandra images of various star cluster types in M51. Above each image the cluster type and B-V color range (from Table [2} is indicated 
along with the number of clusters stacked. Each column corresponds to images stacked in the full (0.3 - 8 keV), hard (2-8 keV), or soft (0.3 - z keV) energy 
bands. The images are 21 pixels per side with each pixel 0.492''. North is up and east is to the left. The source aperture is 6 pixels in diameter while the 
background was calculated outside the larger aperture of diameter 10 pixels. Values in the top left corners of images indicate the 1.3cr upper limit per cluster in 
the stack on the source luminosity in erg s~\ where an asterisk denotes a negative luminosity in the source aperture (below the background level). The grayscale 
denotes the brightest pixel (darkest) relative to the dimmest pixel (lightest). All stacked images across all energy bands correspond to non- detections. 

of the decrease in counts N (fewer clusters stacked), which 
increases the noise (1 / Va^) and hence our upper limits. 

Figure shows the upper limits in erg s"^ plotted against 
the type of cluster that was stacked in the full (blue circles), 
hard (red triangles), and soft (green squares) bands. The re- 
sults from all four methods mentioned above (stacking with 
no modifications as in Figure [6| stacking using the method 
outlined in Figure [7] stacking after removing the 150 bright- 
est (Hof luminosity > 10^^ erg s"^) HII regions from the 
stack, and stacking applying both of the last two methods) 
are presented in Table |4] For all the luminosities calculated, a 
foreground absorption term (assuming a Galactic Wh value of 
1.52 X 10^^ cm"^ and power-law spectrum of F = 1.7) would 
increase our upper limits by ~ 4%. 



E 50 




0.3 0.4 0.5 

Mean of 21 x 21 Region [counts pixel ^] 

Fig. 7. — Histogram representing the mean value of a pixel in the entire 
21x21 region for each young star cluster stacked. Some bins overlap due 
to the large bin size we have chosen for illustration purposes. The long tail 
constitutes star clusters with a large average background value. By removing 
these from the final stacked image we exclude the bright point sources in the 
background region. 



constant (0' - 2', 2' - 3', 3' - 4') and adjusting the source and 
background regions according to PSF size. The result was an 
increase to our upper limits by ~ 15% that is a consequence 



4.1.1. Upper Limits 

An interesting result from Figure [8| is that the young clus- 
ter soft and hard-band limits in panel (a) are above the full- 
band value. Intuitively, the sum of the soft and hard-band 
counts should be equal to the counts in the full-band, but this 
is not true for the luminosity (flux). This discrepancy is a 
consequence of the procedure used to obtain the fluxed im- 
ages in each energy band. From Section 2.1.1| we mentioned 
that a monoenergy was used to calculate the exposure map 
in order to obtain a fluxed image in each energy band. This 
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Fig. 8. — Upper limits in erg plotted against cluster type in the full (blue circles), hard (red triangles), and soft (green squares) bands. Filled shapes correspond 
to positive source luminosities whereas hollow shapes represent negative source luminosities. The values above the x-axes in panels (a) and (c) are the number 
of clusters of each type stacked whereas (b) and (d) show the full, hard, and soft-band numbers from top to bottom respectively. Panel (a) shows the results with 
no modifications (Figure[6j, and (b) shows the results after removing the tail of the histogram (see FigurejTJ from the stacked image, which represented the star 
clusters with the largest average pixel values. Panel (c) corresponds to the removal of the 150 brightest (Ha luminosity > 10^^ erg s~^) HII regions from the stack. 
Panel (d) shows the outcome of applying both methods. The smaller number of clusters present in a stacked image led to larger upper limits since a decrease in 
counts N increases the noise (1 / Va^). 



is not as accurate as using a weighted spectrum file for the 
diff'erent energy bands because the eff'ective area is not con- 
stant throughout the entire energy range. Due to the strong 
energy dependence of eff'ective area, systematic uncertainties 
are introduced that alter the flux value in each pixel from its 
true value. ^ Since the monoenergy used for the hard-band 
was ~ 4.7 keV, all energies up to the 8 keV limit in this 
band would have had their flux underestimated and conse- 
quently introduced a smaller background value. A similar ar- 
gument can be made for the monoenergy of ~ 0.93 keV used 
for the soft-band. These smaller background values led to 
a larger background-subtracted source flux and consequently 
luminosity. The uncertainties in luminosity wer e de termined 
using the conversion factor mentioned in Section [4~T| which is 
the ratio of source luminosity to source counts. Since this fac- 
tor increases with source luminosity it results in much larger 
upper limits. This occurs only for the young clusters because 
we expect there to be more difl'use emission concentrated near 
them. Nevertheless, the variations are small and do not pose 
any significant problems with the analysis. 
The upper limits vary based on numerous characteristics of 

^ |h ttp://c xc. harvard .edu/ciao/threads/spectral.weights/l 
[index . h tml#in tro] 



the energy band, cluster type, and stacking method. In panel 
(c) of Figure [8] the soft-band upper limit for young clusters 
increases above that of the hard-band. This change results 
from the fact we have removed HII regions from the images, 
whic h are known to e mit soft X-rays around kTg = 0.6 ± 0.2 
keV ( Ttillmann et al . [2009) ) . This left mostly harder emission 
features in the soft-band. From our explanation above, this 
should underestimate the background even more and lead to 
an overestimation of the upper limit. 

5. DISCUSSION 
5.1. Upper Limits 

Liu & Li ( 2006]) state that most LMXBs spend their time in 
quiescence at luminosities of < 10^"^ erg s"\ meaning that our 
upper limits are just above that threshold to detect this faint 
population. Our lowest upper limit was found for all clusters 
in the hard band, 6.5 x 10 erg s"^ per cluster. However, this 
is only a l.Scr upper limit and occurs in the hard band, where 
we don't necessarily expect faint XRBs to be found. 

The general trend in Figure [5] follows that the upper limit 
and number of clusters stacked are inversely proportional. All 
our upper limits are below the threshold of 5x10^^ erg s"^ 
which was the minimum X-ray luminosity from the catalogue 
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of point sources (ILllI). Our largest upper limit was from panel 
(a) of Figure [5] for the old cluster population in the full-band, 
9 X 10^^ erg s"\ more than an order of magnitude below the 
threshold luminosity of catalogued point sources. This means 
we have the highest probability of observing bright sources 
(LMXBs) in old clusters compared to other cluster types. Fu- 
ture deeper (750 ks) observations of M51 with Chandra will 
be able to detect faint XRBs in star clusters. 

5.2. Comparisons to the Milky Way 
5.2.1. Globular Clusters 

To put our analysis into context, we need to compare our 
results to those found for other galaxies, such as the Milky 
Way (MW). From a catalogue of LMXBs in the MW (Liu 
^t al. 2007 ) we selected the 13 LMXBs as sociate d with GCs. 
Using distances from the GC catalogue by Harris| ( p^96) [2010 
edition], we find a total X-ray luminosity of GC-LMXBs in 
the MW to be ~ 1.5 x 10 ^^ erg s"^ From the 157 GCs in the 
MW reported by [Harris] ( p^9 6 )[2010 edition] we determine 
the average X-ray luminosity of MW GCs to be ^ 10^^ erg 
s-i. 

From Figures [T] and |4] we assumed that matches 2, 12, 18, 
and 19 are all heavily reddened. In addition, match 20 was as- 
sumed to be an absorbed background source (AGN) because 
it is located far from the extent of the spiral arms and has red 
V - I colors. Therefore only match 4 is likely to be a GC- 
LMXB (even though it is classified as an absorbed source). 
This gives us a total X-ray luminosity of ~ 3.15 x 10^^ erg 
s"\ meaning that the 65 GCs in M51 currently have an aver- 
age X-ray luminosity of ~ 5 x 10^^ erg s"^ This is half the 
value obtained for the average luminosity of MW GCs while 
the total X-ray luminosity is almost an order of magnitude 
smaller. This suggests we have not yet observed a fraction of 
LMXBs in M51. From this approximation we can say that 
either a few (quiescent) bright iq^^"^^ erg s"^) or a large 
population of faint (< 10^^ erg s"^) GC-LMXBs should exist 
in M5 1 . The more likely scenario would be that a combina- 
tion of both (quiescent) bright and faint LMXBs are present in 
M5rs GCs to account for the lower average X-ray luminosity 
compared to the MW. It is also possible that the GC catalogue 
in M51 is incomplete. 

5.2.2. Young Clusters 

The connection betwee n HMXBs and SFR h as been well- 
studied in spiral galaxies (Grim m et al.|2003a|b] ). These stud- 
ies showed that X-ray luminositv increases wit h SFR, mean- 
ing th at the higher SFR of M51 (1.8 M© yr" MColbert et al 
2004| ) would suggest that it contains more XRBs (field and 
cluster) than t he MW (0.68 - 1.45 M p yr'^ [Rob itaille & 
|Whitney"2010). Using the relation from Mine o et al.| ( [2011) , 
we find predicted values of I^xrb (M51) = 4.7 x 10^^ erg s"^ 
and -LxRB (MW) = (1.8 - 3.8) x 10^^ erg s'^ 

However, we want to know the XRBs that are present only 
in clusters. By matching the YCs in the MW with a catalogue 
of HMXBs (within 60 pc), we found the total X-ray lumi- 
nosity of 20 YC-HMXBs to be ~ 2.7 x lO^'^ erg s"^. For the 
MW, the YC population is not comple te past a few kiloparsecs 
(kpc) due to interstellar extinction. Dias et al. j (|2Q02)[2012 
edition] have confirmed ^2100 YCs in the MW, meaning that 
the average X-ray luminosity we would expect from young 
clusters would be ~ 1.3 x 10^"^ erg s"^ . This value is of course 
a crude approximation s ince we observe X-ray point sources 
further away (~ 8 kpc; Grimm 2003 j than we do YCs. Of 



approximately 10 (33%) HMXBs matched to YCs neither the 
HMXB or cluster had distance information and therefore were 
excluded. From Figures [T] and [4] we found numerous X-ray 
point sources that could be classified as HMXBs. They have 
a collective luminosity of ^ 2 x 10^^ erg s"\ which gives an 
average X-ray luminosity of ^ 10^^ erg s"^ for all young clus- 
ters in M5 1 . This value is two orders of magnitude larger than 
that calculated for the MW and likely a result of M5 1 's higher 
SFR. The XR B luminosities calculated using the relation from 
Mineo et al. ([2011) are much larger than those we have cal- 
culated since mey also include field XRBs along with those 
located in clusters. Interestingly, the total X-ray luminosities 
for young and old clusters in M5 1 is approximately identical, 
~ 3.15 X 10^^ erg s'^ for GC-LMXBs and ~ 2.7 x 10^^ erg s'^ 
for YC-HMXBs. This is in spite of the GCs being metal-poor 
and the YCs metal-rich. 

These calculations suff'er from uncertainty due to the tran- 
sient nature of XRBs, meaning that it is possible observations 
were made when some of the brightest XRBs were in a qui- 
escent state. This would skew the number of bright XRBs 
and ultimately produce lower average X-ray luminosities. A 
partial solution to this problem would be to stack multiple ob- 
servations (e.g., as Sivako ff et al. 2008 did for NGC 4697) of 
M5 1 in the hope that any transient behaviour would be aver- 
aged out. Our predictions above for GCs in M5 1 indicate we 
should observe both more (quiescent) bright and faint LMXBs 
coincident with star clusters. This will allow us to perform a 
statistically more meaningful analysis. In addition, the incom- 
pleteness of the YC catalogue in the MW allows for a possibly 
undetected YC-HMXB population. 

5.3. Ejections 

In Section [T] we summarized the methods of ejecting a bi- 
nary from a star cluster. These included natal kicks imparted 
to a NS/BH or fro m binary-binary or three-body interactions. 
From section 3.2 we found that only ~ 1% of M5rs star clus- 
ters have X-ray emission. Besides detecting faint or transient 
X-ray sources coincident with clusters, there are other pro- 
cesses that could account for this small fraction. 
The most favourable among these would be ejections from 



star clusters. Model s and observation s (Rangelo v et al. 
Kaaret et aLl|2004t |Kalogera| |2Q06l JSepinsky et al. 



20TT1 



20051 



showed that binaries can be located > 200 pc from their parent 
cluster. Upon inspection of our stacked images in Figure|6] we 
can see that there are an abundance of point sources spread 
throughout, especially in the stacks of YCs. Our source aper- 
ture of radius 1.5'' covers a range of ~ 120 pc from a cluster 
core. Knowing that most star clusters have cores of only a 
few pc in size, with GCs being slightly larger, an X-ray point 
source located within the cluster should appear directly in the 
centre of an image. It is evident that most of the point sources 
appear outside the background apertures at distances > 200 pc 
from the cluster cores. We cannot rule out the possibility that 
these point sources are XRBs that were in the past located in 
a star cluster. M5rs YCs may have experienced ejections as 
compact objects were created, however, we have no evidence 
for this. 

6. SUMMARY & CONCLUSIONS 

We have presented an X-ray analysis of the star clusters in 
M5 1 to probe the XRB population. We matched 25 star clus- 
ters to 20 X-ray point sources (within 1.5'') and found that the 
majority were representative of HMXBs while those remain- 
ing were characteristic of SNRs. Of the clusters in M51, ~ 
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1% have detected X-ray point sources. Sivakoff et al. (2007 ) 
determined ~ 4% of GCs in 1 1 Virgo Cluster early-type galax- 
ies had X-ray emission, which was larger than the ~ 0.5% we 
found when using the same matching radius of only ^^ The 
discrepancy could have been a result of the short lifetimes of 
HMXBs combined with the larger time period available for 
LMXB formation in GCs. Our stacking of star cluster posi- 
tions in an X-ray image was completed to investigate the low- 
luminosity end of the XLF. Stacking various cluster types in 
the full, hard, and soft energy bands resulted in non-detections 
throughout. The largest upper limit (l.Scr) was found to be 
9.23 X 10^"^ erg s"^ for the average GC in the full-band while 
the lowest upper limit was 6.46 x 10^^ erg s"^ for all clusters in 
the hard-band. The latter limit indicates that any source detec- 
tions would have to consist of faint XRBs. The former limit 
would require that GCs have the largest luminosities in com- 
parison with other cluster types. However, this was biased by 
the small sample size of GCs. 

We implemented a number of methods in order to reduce 
the background level to identify any source signal that was 
present. This included removing bright HIT regions and all 
point sources that were present in the background region. In 
addition, we stacked clusters that were brightest and most 
compact since they are more likely to host XRBs. This pre- 
vented any possible signal from being averaged out by stack- 
ing clusters with no X-ray emission. These methods all re- 
sulted in non-detections. 

The total X-ray luminosities of young and old clusters in 
M51 are 2 x 10^^ erg s"^ and ^ 10^^ erg s"^ respectively. 
Comparing this to clusters in the MW, there should be a pop- 



ulation of undetected transient and faint XRBs in MSl's GC 
population. In contrast, MSl's YC population has more X-ray 
emission than that of the MW's due to its higher SFR. Aver- 
age X-ray luminosities of young and old clusters in M5 1 were 
^ 10^^ erg s"^ and 5 x 10^^ erg s"^ respectively. From in- 
spection of the stacked images it follows that numerous bright 
point sources could have been ejected from clusters in the 
past. If the same is true for faint XRBs then they would rep- 
resent the low-luminosity sources we are searching for. Ulti- 
mately, we conclude that the data is not deep enough to iden- 
tify faint XRBs. With upcoming 750 ks Chandra observa- 
tions of M51 in 2012 September (RL: K. D. Kuntz), we aim 
to classify the low-luminosity XRBs present in star clusters. 
Comparisons with other spiral galaxies and a contrast with el- 
lipticals is necessary to develop a complete understanding of 
faint extragalactic X-ray sources. 
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